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ABSTRACT 

We  predict  the  radio  flux  densities  of  the  extrasolar  planets  in  the  current  census,  making  use  of  an  empirical 
relation — the  radiometric  Bode’s  law — determined  from  the  five  “magnetic”  planets  in  the  solar  system  (the 
Earth  and  the  four  gas  giants).  Radio  emission  from  these  planets  results  from  solar  wind-powered  electron 
currents  depositing  energy  in  the  magnetic  polar  regions.  We  find  that  most  of  the  known  extrasolar  planets 
should  emit  in  the  frequency  range  10-1000  MHz  and,  under  favorable  circumstances,  have  typical  flux  densities 
as  large  as  1  mJy.  We  also  describe  an  initial,  systematic  effort  to  search  for  radio  emission  in  low  radio  frequency 
images  acquired  with  the  Very  Large  Array  (VLA).  The  limits  set  by  the  VLA  images  («300  mJy)  are  consistent 
with,  but  do  not  provide  strong  constraints  on,  the  predictions  of  the  model.  Future  radio  telescopes,  such  as  the 
Low  Frequency  Array  and  the  Square  Kilometer  Array,  should  be  able  to  detect  the  known  extrasolar  planets  or 
place  austere  limits  on  their  radio  emission.  Planets  with  masses  much  lower  than  those  in  the  current  census  will 
probably  radiate  below  10  MHz  and  will  require  a  space-based  array. 

Subject  headings:  magnetic  fields  —  radio  continuum:  stars  —  solar  wind 
Online  material:  machine-readable  table 


1.  INTRODUCTION 

The  past  few  years  have  been  an  exciting  time  as  extrasolar 
planets  have  been  demonstrated  to  be  widespread  and  multiple 
planetary  systems  have  been  found.  The  current  census  now 
numbers  more  than  1 00  extrasolar  planets,  in  over  90  plane¬ 
tary  systems  (Schneider  2003;  Marcy  2003). 

The  vast  majority  of  these  extrasolar  planets  have  been 
detected  via  the  reflex  motion  of  the  host  star.  As  the  existing 
census  shows,  this  method  has  proven  to  be  wildly  successful. 
Nonetheless,  the  reflex  motion  of  the  star  is  a  measure  of  the 
planet’s  gravitational  influence  and  is  necessarily  an  indirect 
detection  of  the  planet.  As  a  consequence,  the  only  property  of 
the  planet  that  one  can  infer  is  its  mass,  and  because  of  the 
mass  function’s  dependence  on  the  inclination  angle  (sin/), 
one  can  infer  only  a  minimum  mass. 

Direct  detection  of  reflected,  absorbed,  or  emitted  radiation 
from  a  planet  allows  for  the  possibility  of  complementary  in¬ 
formation  and  likely  a  more  complete  characterization  of  the 
planet.  The  prototype  of  such  a  direct  detection  is  the  detection 
of  sodium  absorption  lines  in  the  atmosphere  of  the  planet 
orbiting  HD  209458  (Charbonneau  et  al.  2002).  Unfortunately, 
the  incidence  of  transiting  planets  will  always  remain  low 
relative  to  the  total  number  of  planets  known. 

The  Earth  and  gas  giants  of  our  solar  system  are  described 
commonly  as  “magnetic  planets”  because  they  contain  inter¬ 
nal  dynamo  currents  that  are  capable  of  generating  a  planetary- 
scale  magnetic  field.  The  dynamo  currents  themselves  arise 
from  the  rapid  rotation  of  a  conducting  fluid.  The  composi¬ 
tion  of  the  fluid  varies  from  planet  to  planet,  being  liquid  iron 
in  the  Earth,  probably  metallic  hydrogen  in  Jupiter  and  Saturn, 


and  probably  a  salty  ocean  in  Uranus  and  Neptune.  In  turn, 
these  planetary  magnetic  fields  are  immersed  in  the  high-speed 
electrical  gas  emitted  by  the  Sun.  Via  a  coupling  between  the 
solar  wind  and  the  planetary  magnetic  field,  all  of  these  mag¬ 
netic  planets  produce  radio  emission. 

Some  of  the  extrasolar  planets  are  also  magnetic  planets. 
Shkolnik  et  al.  (2003)  have  detected  a  modulation  in  the  Ca  n 
H  and  K  lines  of  HD  179949  with  a  periodicity  that  is  that  of 
the  planetary  orbit.  They  interpret  this  as  a  magnetic  interac¬ 
tion  between  the  star  and  planet,  although  there  is  no  con¬ 
straint  as  yet  on  the  magnetic  field  strength  of  the  planet. 
Other  stars  in  their  sample  also  appear  to  have  modulations  in 
the  Ca  ii  H  and  K  lines,  although  it  is  not  yet  clear  that  these 
modulations  are  periodic  with  the  orbital  period  of  their 
planets,  as  is  the  case  for  HD  179949. 

By  analogy  to  solar  system  planets,  there  have  been  a 
number  of  suggestions  and  searches  for  natural  radio  emission 
from  extrasolar  planets.  Early  searches  were  conducted  by 
Yantis  et  al.  (1977)  and  Lecacheux  (1991).  Because  of  the  lack 
of  any  known  extrasolar  planets  at  the  time,  these  were  nec¬ 
essarily  blind  searches.  More  recently,  Farrell  et  al.  (1999)  and 
Zarka  et  al.  (2001),  building  on  empirical  relations  derived 
for  solar  system  planets  and  calibrated  by  spacecraft  fly-bys, 
have  made  specific  predictions  for  the  radio  emission  from  the 
known  extrasolar  planets. 

Detection  of  radio  emission  from  an  extrasolar  planet  would 
constitute  a  direct  detection  and  can  yield  fundamental  in¬ 
formation  about  the  planet.  First,  a  measurement  of  the  radio 
emission  is  directly  indicative  of  the  polar  magnetic  field 
strength  at  the  planet.  For  example,  the  high-frequency  cutoff 
of  Jovian  decametric  bursts  (~40  MHz)  is  interpreted  as  being 
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due  to  the  Jovian  polar  magnetic  field  strength,  which  allowed 
an  estimate  for  the  strength  of  the  Jovian  magnetic  field  nearly 
20  yr  prior  to  the  first  in  situ  magnetic  field  observations. 
In  turn,  of  course,  the  presence  of  a  magnetic  field  provides 
a  rough  measure  of  the  composition  of  the  planet,  insofar  as 
it  requires  the  planet’s  interior  to  have  a  conducting  fluid. 
Combined  with  an  estimate  of  the  planet’s  mass,  one  could 
deduce  the  composition  of  the  fluid  by  analogy  to  the  solar 
system  planets  (liquid  iron  vs.  metallic  hydrogen  vs.  salty 
ocean). 

Second,  the  periodic  nature  of  the  radio  emission  has  been 
used  to  define  precisely  the  planetary  rotation  periods  of  all 
of  the  gas  giant  planets  in  the  solar  system.  As  the  magnetic 
field  is  presumed  to  be  tied  to  the  interior  of  the  planet,  it 
provides  a  more  accurate  measure  of  the  planet’s  rotation  rate 
than  atmospheric  phenomena  such  as  clouds.  For  instance, 
the  rotation  period  of  Neptune  was  determined  initially  by 
observations  of  the  differentially  rotating  cloud  tops  but  then 
was  redefined  after  detection  of  the  Neptunian  radio  emission 
(Lecacheux  et  al.  1993). 

Finally,  testing  the  extent  to  which  solar  system  models  of 
magnetic  fields  can  be  applied  to  extrasolar  planets  may  have 
important  implications  for  assessing  the  long-term  “habit¬ 
ability”  of  any  terrestrial  planets  found  in  the  future.  The  im¬ 
portance  of  a  magnetic  field  is  that  it  deflects  incident  cosmic 
rays.  If  these  particles  reach  the  surface  of  an  otherwise  hab¬ 
itable  planet,  they  may  cause  severe  cellular  damage  and  dis¬ 
ruption  of  genetic  material  to  any  life  on  its  surface  or  may 
prevent  life  from  arising  at  all.  A  secondary  importance  of  a 
magnetic  field  is  that  it  can  prevent  the  planet’s  atmosphere 
from  being  eroded  by  the  stellar  wind  (Mitchell  et  al.  2001); 
this  process  is  thought  to  be  a  contributing  factor  to  the  relative 
thinness  of  the  Martian  atmosphere.  (This  is  of  course  unlikely 
to  be  an  issue  for  extrasolar  giant  planets.) 

This  paper  extends  and  updates  initial  attempts  to  predict  the 
radio  emission  of  extrasolar  giant  planets  by  applying  the  solar 
system-determined  “radiometric  Bode’s  law”  to  the  current 
census  of  extrasolar  planets.  In  §  2  we  summarize  briefly  the 
various  forms  of  the  radiometric  Bode’s  law  and  the  form  we 
have  adopted  for  use,  in  §  3  we  summarize  the  current  state  of 
the  extrasolar  planet  census  and  predict  the  “burst”  emission 
levels  and  frequencies  for  radio  emission  from  these  planets,  in 
§  4  we  use  the  initial  observations  from  a  74  MHz  sky  survey 
to  confront  the  model  predictions  with  data,  and  in  §  5  we 
assess  how  existing  and  future  radio  telescopes  may  be  able  to 
detect  radio  emission  from  extrasolar  planets  and  present  our 
conclusions. 

2.  THE  RADIOMETRIC  BODE’S  LAW 

All  of  the  planets  in  the  solar  system  with  substantial 
magnetic  fields  are  strong  emitters  of  coherent  cyclotron 
emission  at  radio  wavelengths.  The  emission  itself  arises  from 
energetic  (keV)  electrons  propagating  along  magnetic  field 
lines  into  active  auroral  regions  (Gumett  1974;  Wu  &  Lee 
1979).  The  source  of  the  electron  acceleration  to  high  energies 
is  ultimately  a  coupling  between  the  incident  solar  wind  and 
the  planet’s  magnetic  field.  The  precise  details  are  not  well 
understood,  but  the  coupling  arises  from  magnetic  field  re¬ 
connection  in  which  the  magnetic  field  embedded  in  the  solar 
wind  and  the  planetary  magnetic  field  cancel  at  their  interface, 
thereby  energizing  the  plasma.  Energetic  electrons  in  the  en¬ 
ergized  plasma  form  a  current  flow  planetward  along  the 
planet’s  magnetic  field  lines,  with  the  lines  acting  effectively 


like  low-resistance  wires.  The  energy  in  these  magnetic  field- 
aligned  electric  currents  is  deposited  in  the  upper  polar  at¬ 
mosphere  and  is  responsible  for  the  visible  aurora.  Besides 
auroral  emissions  at  visible  wavelengths  (e.g.,  northern  lights), 
about  10~5  to  1CU6  of  the  solar  wind  input  power  is  converted 
to  escaping  cyclotron  radio  emission  (Gumett  1974).  The  au¬ 
roral  radio  power  from  the  Earth  is  of  order  10s  W,  while  for 
Jupiter  it  is  of  order  10105  W. 

The  cyclotron  emission  process  described  is  generic,  and 
it  is  typically  most  efficient  in  the  auroral  region,  about  1-3 
planetary  radii  in  altitude.  Specific  details  of  the  cyclotron 
emission  process  do  differ  from  planet  to  planet,  depending 
on  secondary  effects  such  as  the  polar  source  location  and 
beaming  angles,  which  in  turn  depend  on  the  planet’s  magnetic 
field  topology  and  low-altitude  atmospheric  density.  None¬ 
theless,  applicable  to  all  of  the  magnetic  planets  is  a  mac¬ 
roscopic  relationship  called  the  “radiometric  Bode’s  law” 
relating  the  incident  solar  wind  power,  the  planet’s  magnetic 
field  strength,  and  the  median  emitted  radio  power.  Various 
investigators  (Desch  &  Barrow  1984;  Desch  &  Kaiser  1984; 
Desch  &  Rucker  1985;  Barrow  et  al.  1986;  Rucker  1987; 
Desch  1988;  Millon  &  Goertz  1988)  have  found 

^rad  =  <v,  (1) 

where  Pad  is  the  median  emitted  radio  power,  Psw  is  the  in¬ 
cident  solar  wind  power,  e  ~  10-5  to  10-6  is  the  efficiency  at 
which  the  solar  wind  power  is  converted  to  emitted  radio 
power,  and  x  ~  1 .  The  values  of  the  coefficient  e  and  the 
power-law  index  x  measured  by  various  investigators  have 
differed  over  the  years  on  the  basis  of  the  quantity  of  data 
available;  the  early  determinations  were  based  on  fly-bys  of 
only  Jupiter  and  Saturn  and  measurements  of  the  Earth’s  radio 
emission,  while  later  determinations  include  data  from  the 
Voyager  2  Uranus  and  Neptune  fly-bys  as  well.  In  general,  the 
more  recent  work  indicates  higher  efficiencies  than  the  earlier 
work  (Zarka  et  al.  1997). 

Farrell  et  al.  (1999)  and  Zarka  et  al.  (2001)  extrapolated 
this  relation  to  extrasolar  planetary  systems  by  making  use  of 
two  scaling  relations  and  adopting  more  recent  estimates  for  e 
and  x.  The  incident  solar  wind  power  depends  on  the  ram 
pressure  of  the  solar  wind  and  the  cross-sectional  area  pre¬ 
sented  by  the  magnetosphere  to  the  solar  wind.  The  ram 
pressure  of  the  solar  wind  depends  on  the  star-planet  distance 
d,  while  the  size  of  the  planet’s  magnetosphere  is  related 
to  its  magnetic  moment.  Blackett  (1947)  demonstrated  that 
planetary-scale,  dynamo-driven  magnetic  fields  are  scalable, 
with  a  functional  dependence  on  body  mass,  radius,  and  an¬ 
gular  frequency;  and  the  presence  of  planetary-scale  magnetic 
fields  around  the  gas  giants  shows  that  his  analysis  does  not 
depend  on  the  actual  composition  of  the  conducting  fluid  in 
the  planet’s  interior.  Blackett’s  law  relates  a  planet’s  magnetic 
moment  (its  surface  field  times  its  radius  cubed,  BR3)  to  its 
rotation  rate  and  mass  by 

(2) 

where  u>  is  the  planet’s  rotation  rate  and  M  is  its  mass.  Since 
that  time  and  with  the  inclusion  of  the  Voyager  magnetic  field 
observations  at  the  outer  planets,  numerous  planetary  mag¬ 
netic  field  scaling  laws  have  been  derived,  each  expanding  on 
but  maintaining  some  similarity  to  Blackett’s  law. 
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Combining  these  various  scaling  relations,  Farrell  et  al. 
(1999)  and  Zarka  et  al.  (2001)  found  that  the  planet’s  median 
radio  power  driven  by  the  stellar  wind  is 
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We  have  normalized  all  quantities  to  those  of  Jupiter.  Farrell 
et  al.  (1999,  Appendix  A)  discussed  slight  differences  to  this 
radiometric  Bode’s  law;  the  differences  result  from  the  sta¬ 
tistical  spread  in  the  various  (solar  system)  planets’  magnetic 
moments  and  amount  to  slightly  different  exponents  and/or  a 
different  coefficient. 

We  stress  various  aspects  of  this  radiometric  Bode’s  law. 
First,  it  is  grounded  in  in  situ  measurements  from  spacecraft 
fly-bys  of  the  gas  giants,  as  well  as  measurements  of  the  Earth’s 
radio  emission.  In  particular,  Farrell  et  al.  (1999)  considered 
two  forms  of  the  radiometric  Bode’s  law,  one  incorporating 
measurements  of  only  Jupiter  and  Saturn  and  the  other  incor¬ 
porating  all  five  magnetic  planets.  The  former  gives  lower 
power  levels  (by  a  factor  of  100);  in  equation  (3)  we  have 
adopted  the  latter  because  it  is  based  on  a  complete  census  of 
magnetic  planets  in  the  solar  system. 

Second  is  the  importance  of  the  planet-star  distance.  As 
Zarka  et  al.  (2001)  show,  the  radio  power  from  the  Earth  is 
larger  than  that  from  Uranus  or  Neptune,  even  though  both  of 
those  planets  have  magnetic  moments  approximately  50  times 
larger  than  that  of  the  Earth. 

Third,  the  power  levels  estimated  by  equation  (3)  are  for 
emission  into  a  solid  angle  of  47r.  Emission  solid  angles  for 
solar  system  planets  span  a  large  range.  Considering  Jovian 
decametric  emission  at  macroscopic  scales,  emission  occur¬ 
rence  is  as  Cl  ~  7 r.  In  estimating  the  flux  density,  the  relevant 
quantity  is  P/Cl.  (See  eq.  [5]  below  and  the  discussion  in  Farrell 
et  al.  1999.)  Assuming  a  smaller  emission  solid  angle  would 
mean  reducing  the  nominal  power  levels  of  equation  (3)  by 
the  factor  Cl/Air.  (This  factor  also  explains  the  difference  be¬ 
tween  the  nominal  power  level  for  Jupiter  predicted  by  eq.  [3] 
and  the  level  quoted  in  §  1 .) 

Fourth,  the  radiometric  Bode’s  law  of  equation  (3)  describes 
the  median  emitted  power  from  the  magnetospheres  of  the 
Earth  and  all  of  the  solar  gas  giants,  including  the  non-Io- 
driven  Jovian  decametric  radio  emission  (Zarka  et  al.  2001). 
Planetary  magnetospheres  tend  to  act  as  “amplifiers”  of  the 
incident  solar  wind,  so  that  an  increase  in  the  solar  wind  ve¬ 
locity  (and  therefore  incident  pressure)  leads  to  geometrically 
higher  emission  levels.  This  effect  is  observed  at  all  of  the 
magnetized  planets.  For  instance,  the  Saturnian  radio  flux  has 
been  observed  to  increase  by  factors  of  20  because  of  changes 
in  the  solar  wind  velocity  (Desch  &  Rucker  1985).  At  the 
Earth,  the  radio  power  has  been  observed  to  increase  by  a 
factor  of  10  for  every  200  km  s_1  increase  in  the  solar  wind 
velocity  (Gallagher  &  D’Angelo  1981,  Fig.  5);  the  observed 
range  of  solar  wind  velocities  at  the  Earth  is  at  least  400  km 
s_1  (from  400  to  800  km  s_1).  The  radio  power  levels  from  a 
planet  can  exceed  that  predicted  by  equation  (3)  by  factors  of 
100  and  possibly  more. 

Finally,  an  implicit  assumption  in  our  analysis  is  that  the 
mass-loss  rates  for  stars  hosting  planets  are  similar  to  that  of 
the  Sun.  The  stellar  wind  pressure  incident  on  a  planetary 
magnetosphere  is  Psw  c x  pF3,  where  p  is  the  density  of  the 
stellar  wind  and  V  is  its  velocity.  Stars  with  higher  (lower) 
mass-loss  rates  than  the  Sun  will  have  higher  (lower)  stellar 


wind  densities  and  should  drive  higher  (lower)  planetary  radio 
emissions.  Differences  in  the  mass-loss  rate  from  star  to  star 
could  lead  to  statistical  variations  in  the  power  levels  we 
predict,  although  these  should  be  relatively  less  important  than 
variability  in  stellar  wind  velocities.  (Work  is  in  progress  to 
utilize  what  is  known  about  the  various  stellar  ages  and  ac¬ 
tivities  to  relax  this  assumption.) 

In  order  for  gyrating  electrons  to  generate  radio  emission 
that  escapes  the  auroral  region,  the  Doppler-shifted  electron 
cyclotron  frequency  must  couple  directly  to  a  magnetoplasma 
mode.  This  coupling  can  only  occur  in  low-density,  high 
magnetic  field  conditions  where  the  local  electron  plasma-to- 
cyclotron  frequency  ratio  is  much  less  than  1.  The  regions  just 
above  the  polar  ionosphere  satisfy  this  coupling  criteria. 
Hence,  the  characteristic  emission  frequency  can  be  related  to 
the  planetary  magnetic  dipole  moment  (eq.  [2])  by 

23-s  mK£)  (£fs/  <4» 


In  order  to  be  detected  from  ground-based  telescopes,  vc  must 
exceed  the  terrestrial  ionospheric  cutoff  (defined  by  the  peak 
ionospheric  plasma  frequency),  which  is  typically  3-10  MHz. 
Farrell  et  al.  (1999)  showed  that  vc  exceeded  the  Earth’s 
ionospheric  cutoff  (~3  MHz)  for  some  planets  known  at  the 
time  of  the  publication  of  their  paper,  with  the  planet  around 
HD  114762  having  vc  «  1000  MHz. 

Just  as  the  statistical  variations  in  the  solar  system  planetary 
magnetic  moments  contribute  to  variations  in  the  radiometric 
Bode’s  law,  there  is  also  a  statistical  spread  in  the  exponents 
and/or  coefficients  for  the  emission  frequency.  Compared  to 
the  radiated  power  predicted  by  equation  (3),  the  statistical 
variations  in  the  predicted  emission  frequency  are  relatively 
more  important.  The  radiated  power  levels  can  be  exponen¬ 
tially  modified  by  the  solar  wind  loading,  whereas  the  emis¬ 
sion  frequency  is  not.  In  order  to  account  for  the  statistical 
variations  in  estimates  of  the  magnetic  moment  derived  from 
Blackett’s  law,  we  assume  that  the  actual  emission  frequency 
could  be  larger  or  smaller  by  a  factor  of  3  than  the  nominal 
value  of  equation  (4).  (See  also  Appendix  A  of  Farrell  et  al. 
1999.) 

The  predicted  flux  density  of  an  extrasolar  planet  is  then 


P 

AuClD2 ' 


(5) 


Farrell  et  al.  (1999)  applied  equations  (2) — (5)  to  the  nine  then- 
known  extrasolar  planets.  The  distances  to  the  host  stars  are 
known.  The  planetary  mass  and  distance  to  its  parent  star  are 
measured  or  constrained  from  the  optical  observations.  As¬ 
suming  an  emission  cone  of  4-7T  sr,  A v  =  vc/2,  and  a  rotation 
rate  comparable  to  that  of  Jupiter,  they  derived  the  resulting 
flux  densities. 

We  do  the  same  but  for  the  currently  much  larger  sample. 

3.  THE  EXTRASOLAR  PLANET  CENSUS 

Since  the  initial  work  by  Farrell  et  al.  (1999),  the  extrasolar 
planetary  census  has  grown  by  over  an  order  of  magnitude. 
We  used  the  census  of  extrasolar  planets  contained  in  “The 
Extrasolar  Planets  Encyclopaedia”  (Schneider  2003).  We  used 
only  the  set  of  reasonably  secure  detections  of  planets  around 
main-sequence  stars  as  of  2003  July  1.  The  total  number  of 
planets  examined  was  118  in  102  planetary  systems. 


514 


LAZIO  ET  AL. 


Vol.  612 


TABLE  1  TABLE  1 — Continued 

Exoplanet  Radio  Emission 


Star 

ui/u>  J 

logAad 

(W) 

A* 

(GRj) 

Vc 

(MHz) 

log  5 

(Jy) 

OGLE-TR-56 . 

0.3 

16.95 

2.6 

14.8 

-4.3 

HD  73256  . 

0.2 

16.77 

3.7 

20.8 

-1.4 

HD  83443  . 

0.1 

15.79 

0.3 

1.4 

-1.4 

HD  46375  . 

0.1 

15.48 

0.1 

0.6 

-1.1 

HD  179949 . 

0.1 

16.28 

1.1 

5.9 

-1.1 

HD  187123  b . 

0.1 

15.88 

0.4 

2.1 

-1.6 

HD  209458 . 

0.1 

15.95 

0.5 

2.9 

-1.6 

HD  75289  . 

0.1 

15.65 

0.2 

1.3 

-1.2 

BD  -10°  3166 . 

0.1 

15.73 

0.3 

1.6 

HD  76700 . 

0.1 

15.13 

0.1 

0.3 

-1.7 

T  Boo  . 

0.1 

16.93 

10.7 

59.9 

-1.0 

51  Peg . 

0.1 

15.57 

0.2 

1.2 

-0.6 

HD  49674 . 

0.1 

14.66 

0.0 

0.1 

-1.4 

v  And  b . 

0.09 

15.67 

0.4 

2.2 

-0.7 

v  And  c . 

1.0 

15.25 

12.1 

68.0 

-2.6 

v  And  d . 

1.0 

14.87 

38.0 

213 

-3.5 

HD  168746 . 

0.07 

14.86 

0.0 

0.3 

-1.6 

HD  217107  b . 

0.06 

15.76 

0.7 

4.0 

-1.8 

HD  68988  . 

0.07 

16.02 

1.6 

8.8 

-2.2 

HD  162020 . 

0.05 

17.06 

31.7 

178 

-2.0 

HD  130322 . 

0.04 

15.36 

0.4 

2.0 

-1.7 

HD  108147 . 

0.04 

14.67 

0.1 

0.4 

-1.8 

G1  86 . 

1 

17.08 

42.3 

237 

-1.2 

55  Cnc  b . 

1 

16.15 

3.1 

17.6 

-1.1 

55  Cnc  c . 

1 

14.84 

0.3 

1.7 

-1.4 

55  Cnc  d . 

1 

14.32 

43.2 

242 

-4.1 

HD  38529  b . 

1 

16.03 

2.8 

15.5 

-2.2 

HD  38529  c . 

1 

15.31 

290.4 

1626 

-4.9 

HD  195019 . 

1 

16.83 

32.8 

184 

-1.8 

HD  6434 . 

1 

15.64 

1.2 

6.9 

-2.2 

HD  192263  . 

1 

15.88 

2.4 

13.6 

-1.6 

G1  876  b . 

1 

16.20 

12.1 

68.0 

-0.8 

G1  876  c . 

1 

15.83 

1.6 

8.9 

-0.2 

p  CrH . 

1 

15.82 

4.9 

27.6 

-1.8 

HD  74156  b . 

1 

15.90 

8.8 

49.4 

-3.2 

HD  74156  c . 

1 

15.05 

120.7 

676 

-5.2 

HD  3651  . 

1 

14.69 

0.3 

1.6 

-1.4 

HD  168443  b . 

1 

16.75 

112.8 

631 

-2.9 

HD  168443  c . 

1 

15.66 

476.7 

2670 

-4.6 

HD  121504 . 

1 

15.47 

3.5 

19.4 

-2.9 

HD  1 789 1 1 B . 

1 

16.60 

90.1 

504.4 

-3.2 

HD  16141 . 

1 

14.59 

0.3 

1.8 

-2.6 

70  Vir . 

1 

16.42 

97.5 

546 

-2.8 

HD  80606 . 

1 

16.03 

32.4 

182 

-3.5 

HD  219542B . 

1 

14.59 

0.6 

3.2 

-3.2 

HD  216770 . 

1 

15.08 

2.3 

13.0 

-3.0 

GJ  3021 . 

1 

15.94 

31.0 

174 

-2.6 

HD  52265  . 

1 

15.31 

5.1 

28.8 

-2.8 

HD  37124  b . 

1 

15.01 

2.6 

14.6 

-3.0 

HD  37124  c . 

1 

14.22 

5.7 

31.9 

-4.1 

HD  73526  . 

1 

15.67 

26.2 

147 

-4.3 

HD  8574 . 

1 

15.40 

16.0 

89.5 

-3.6 

HD  104985  . 

1 

15.98 

90.3 

505 

-4.5 

HD  134987  . 

1 

15.18 

9.0 

50.4 

-3.1 

HD  169830  b . 

1 

15.51 

23.5 

132 

-3.5 

HD  169830  c . 

1 

14.53 

17.2 

96.3 

-4.3 

HD  40979 . 

1 

15.59 

31.0 

174 

-3.5 

HD  150706 . 

1 

14.89 

4.2 

23.5 

-3.1 

HD  12661  b . 

1 

15.36 

16.8 

94.3 

-3.5 

HD  12661  c . 

1 

14.35 

8.9 

49.9 

-4.3 

HD  89744 . 

1 

15.98 

112.8 

631 

-3.8 

HR  810 . 

1 

15.27 

16.2 

90.9 

-2.8 

HD  92788  . 

1 

15.56 

39.5 

221 

-3.6 

HD  142 . 

1 

14.76 

4.2 

23.5 

-3.0 

Star 

uj/uj] 

logAad 

(W) 

(G  R,) 

Vc 

(MHz) 

log  5 

Uy) 

HD  177830  . 

1 

14.89 

6.3 

35.5 

-4.0 

HD  128311 . 

1 

15.30 

21.0 

118 

-3.0 

HD  28185  . 

1 

15.73 

76.4 

428 

-3.9 

HD  108874 . 

1 

14.99 

9.7 

54.2 

-4.2 

HD  142415  . 

1 

15.02 

10.5 

58.6 

-3.6 

HD  4203 . 

1 

14.97 

9.6 

53.6 

-4.3 

HD  210277 . 

1 

14.81 

6.0 

33.7 

-3.2 

HD  82943  b . 

1 

14.93 

9.5 

53.1 

-3.5 

HD  82943  c . 

1 

14.89 

3.4 

19.0 

-3.0 

HD  27442  . 

1 

14.84 

7.6 

42.7 

-3.1 

HD  114783 . 

1 

14.61 

4.1 

23.1 

-3.2 

HD  20367 . 

1 

14.63 

4.7 

26.3 

-3.4 

HD  147513  . 

1 

14.59 

4.2 

23.5 

-2.8 

HD  19994 . 

1 

14.97 

13.3 

74.7 

-3.4 

HD  65216 . 

1 

14.73 

6.8 

37.8 

-3.7 

HD  4 1004 A . 

1 

15.04 

16.8 

94.3 

-4.0 

HIP  75458 . 

1 

15.79 

152.8 

856 

-3.9 

HD  222582 . 

1 

15.48 

63.7 

357 

-4.1 

HD  160691  . 

1 

14.78 

10.2 

57.0 

-3.1 

HD  23079  . 

1 

15.01 

19.9 

111 

-3.9 

HD  141937 . 

1 

15.77 

185.3 

1038 

-4.1 

HD  47536 . 

1 

15.34 

60.6 

339 

-5.1 

HD  114386 . 

1 

14.41 

4.1 

23.1 

-3.6 

HD  4208  . 

1 

14.26 

3.0 

16.6 

-3.8 

16  Cyg  B . 

1 

14.61 

8.3 

46.2 

-3.5 

HD  10697 . 

1 

15.36 

97.3 

545 

-4.1 

7  Cep . 

1 

14.52 

9.1 

50.9 

-3.1 

HD  213240 . 

1 

15.12 

51.5 

289 

-4.3 

HD  111232 . 

1 

15.43 

128.8 

722 

-4.1 

HD  1 14729 . 

1 

14.12 

3.0 

16.9 

-4.0 

47  UMa  b . 

1 

14.77 

19.9 

111 

-3.3 

47  UMa  c . 

1 

13.67 

2.7 

14.9 

-3.5 

HD  10647 . 

1 

14.32 

5.5 

30.6 

-3.4 

HD  2039 . 

1 

15.14 

63.5 

355 

-5.1 

HD  190228 . 

1 

15.09 

61.2 

343 

-4.9 

HD  50554  . 

1 

15.06 

59.4 

333 

-4.2 

HD  136118 . 

1 

15.57 

260.9 

1461 

-4.8 

HD  196050 . 

1 

14.75 

26.2 

147 

-4.5 

HD  30177 . 

1 

15.26 

126.1 

706 

-4.8 

HD  106252 . 

1 

15.19 

102.8 

576 

-4.5 

HD  216435  . 

1 

14.29 

8.2 

45.7 

-4.2 

HD  216437 . 

1 

14.49 

14.5 

81.0 

-4.0 

HD  23596  . 

1 

15.19 

112.5 

630 

-4.8 

14  Her . 

1 

14.94 

59.4 

333 

-3.9 

HD  72659 . 

1 

14.47 

20.0 

112 

-4.8 

e  Eri . 

1 

13.83 

3.3 

18.3 

-2.2 

HD  39091  . 

1 

15.26 

204.8 

1147 

-4.2 

HD  33636  . 

1 

15.39 

183.4 

1027 

-4.3 

Gl  777A . 

1 

13.82 

6.8 

37.8 

-3.9 

Notes. — Planets  are  arranged  roughly  in  order  of  increasing  semimajor 
axis.  The  predicted  flux  densities  are  the  burst  values.  Table  1  is  also  available 
in  machine-readable  form  in  the  electronic  edition  of  the  Astrophysical 
Journal. 

The  required  quantities  for  applying  equation  (3)  are  the 
planet’s  mass,  distance  from  its  primary,  and  rotation  rate.  At 
the  time  we  extracted  our  list  of  planets  from  “The  Extrasolar 
Planets  Encyclopaedia,”  almost  all  of  the  known  planets  had 
been  found  via  the  radial  velocity  method.  This  method  yields 
a  lower  limit  to  the  planet’s  mass  (i.e.,  M  sin  i). 

Strictly,  with  the  exception  of  the  planets  orbiting  OGLE- 
TR-56  and  HD  209458,  for  which  transits  have  been  detected, 
the  flux  density  levels  that  we  calculate  are  lower  limits.  If  the 
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Emission  Frequency  (MHz) 

Fig.  1. — Predicted  “burst”  flux  densities  for  106  known  extrasolar  planets 
vs.  the  characteristic  emission  frequency  based  on  the  radiometric  Bode’s  law 
and  Blackett’s  law(  eqs.  [3] — [5]).  The  horizontal  bars  indicate  the  assumed 
ranges  for  the  emission  frequencies,  allowing  the  statistical  variations  from 
Blackett’s  law  in  the  estimated  planetary  magnetic  moments.  The  expected 
burst  flux  densities  are  obtained  by  assuming  that  increases  of  roughly  a  factor 
of  100  can  be  obtained  by  larger  stellar  wind  loading  of  the  planet’s  magne¬ 
tosphere.  The  vertical  dashed  line  indicates  the  approximate  cutoff  frequency 
for  Jupiter. 

emitted  radiation  is  nearly  isotropic  and  assuming  that  the 
inclinations  of  the  orbits  are  random,  the  typical  error  intro¬ 
duced  is  only  a  factor  of  2. 

There  have  been  attempts  to  set  upper  limits  on  planetary 
masses,  e.g.,  using  tidal  constraints  (Trilling  2000).  We  have 
not  made  use  of  these  upper  limits  because  in  many  cases  they 
introduce  additional  uncertainties  based  on  the  characteristics 
of  the  host  star  (e.g.,  its  age). 

The  radial  velocity  method  also  determines  the  semimajor 
axis  and  eccentricity  of  the  planet’s  orbit.  We  use  the  semimajor 
axis  as  the  distance  d  in  equation  (3).  For  those  planets  in  ec¬ 
centric  orbits,  we  might  expect  the  radio  emission  to  be  peri¬ 
odic  with  the  planet’s  orbital  period,  increasing  near  periastron 
and  decreasing  near  apastron.  The  most  extreme  eccentricities 
in  our  sample  are  approximately  0.7  (HIP  75458,  HD  222582, 
and  HD  2039).  We  therefore  expect  a  range  of  possible  flux 
densities  of  possibly  20.  For  a  planet  with  a  more  typical 
eccentricity,  e  «  0.1,  the  range  is  only  a  factor  of  1.5  or  so. 

We  have  no  information  on  the  rotation  rates  of  these 
planets.  We  have  therefore  assumed  that  those  planets  with 
semimajor  axes  larger  than  0. 1  AU  have  rotation  rates  equal  to 
that  of  Jupiter  (10  hr),  while  those  with  semimajor  axes  less 
than  0.1  AU  are  tidally  locked,  with  rotation  rates  equal  to 
their  orbital  periods  (Goldreich  &  Soter  1966;  Trilling  2000). 

For  tidally  locked  planets  we  have  also  assumed  a  radius  of 
1.25Rj  in  calculating  the  characteristic  frequency.  This  value  is 
both  predicted  based  on  the  stellar  insolation  of  such  planets 
(Guillot  et  al.  1996)  and  is  close  to  the  observed  value  for  the 
planet  orbiting  HD  209458  (Charbonneau  et  al.  2002). 

Table  1  presents  the  full  list  of  planets  in  our  sample  along 
with  the  expected  emission  frequencies  and  flux  density  levels. 
Figure  1  presents  the  expected  flux  densities  versus  the  emis¬ 
sion  frequency  in  a  graphical  form .  Anticipating  our  analysis  of 
§  4,  we  follow  the  convention  of  Farrell  et  al.  (1999)  and  in¬ 
crease  the  predicted  flux  density  levels  by  a  factor  of  1 00  in  order 
to  account  for  variability  driven  by  the  incident  stellar  wind. 

The  trend  in  Figure  1  of  increasing  emission  frequency  and 
decreasing  flux  density  is  real  and  reflects  two  effects.  First, 


the  lower  envelope  reflects  a  selection  effect.  A  low  flux 
density  and  small  emission  frequency  results  from  a  low- 
mass  planet  in  a  large  orbit.  These  planets  cannot  be  detected 
with  the  current  detection  methodology.  Second,  the  upper 
envelope  reflects  the  well-known  deficit  of  planets  with  both 
large  masses  and  small  semimajor  axes.  Even  if  high-mass 
planets  with  close  orbits  did  exist,  however,  they  probably 
would  be  tidally  locked  (as  we  have  assumed  here).  The 
rotation  rate  also  determines  the  strength  of  the  magnetic 
field  (eq.  [4]),  so  tidally  locked  planets  probably  do  not  ra¬ 
diate  at  high  emission  frequencies.  However,  Zarka  et  al. 
(2003)  have  suggested  that  such  “hot  Jupiters”  may  radiate 
by  the  conversion  of  stellar  magnetic  pressure  (which  is  large 
in  close  to  the  star)  into  electromotive  forces,  currents  and 
radio  emission. 

We  have  assumed  that  the  objects  discussed  here  (Table  1) 
are  in  fact  planets.  Using  Hipparcos  astrometric  data,  Han  et  al. 
(2001)  argue  that  the  inclinations  of  the  orbits  of  these  objects 
imply  that  sin  i  is  small.  A  small  sin  i  would  imply  that  the 
derived  masses  of  these  objects  would  be  larger,  so  that  their 
magnetic  moments  would  be  larger,  but  also  that  their  dis¬ 
tances  from  their  stars  would  be  larger,  so  that  the  stellar  wind 
loading  of  the  magnetosphere  would  be  less.  The  similar 
dependences  for  M  and  d  in  equation  (3)  mean  that  this  impact 
of  small  sin  i  would  nearly  cancel.  However,  if  sin  i  is  small 
enough,  as  Han  et  al.  (2001)  suggest  for  a  number  of  objects, 
they  might  no  longer  be  classified  as  planets  but  as  brown 
dwarfs.  Berger  et  al.  (2001)  detected  radio  emission  from  one 
brown  dwarf,  indicating  that  it  is  likely  to  be  a  magnetic  body 
as  well.  If  so,  the  radiometric  Bode’s  law  and  our  analysis 
might  still  apply  to  them. 

We  now  discuss  selected  objects  in  greater  detail. 

OGLE-TR-56. — The  planet  orbiting  this  star  was  found  from 
its  transits  of  its  parent  star  (Konacki  et  al.  2003).  In  many 
respects  it  is  not  much  different  than  many  of  the  other  “close- 
in”  planets.  Its  flux  density  is  roughly  3  orders  of  magnitude 
lower,  though,  because  of  its  much  greater  distance  (1500  vs.  a 
more  typical  10  pc).  This  result  is  likely  to  be  a  general  result 
as  future  planetary  searches  move  deeper  into  the  Galaxy:  the 
detection  of  extrasolar  planetary  radio  emission  is  most  likely 
to  come  from  planets  orbiting  nearby  stars. 

HD  179949. — A  possible  magnetic  interaction  between 
the  planet  and  the  star  has  been  detected  for  this  star  (Shkolnik 
et  al.  2003). 

r  Boo,  HD  209458. — These  stars  show  Ca  ii  H  and  K  ac¬ 
tivity  similar  to  that  of  HD  179949.  However,  a  clear  link 
between  the  planetary  orbital  period  and  the  modulation  of  the 
Ca  ii  H  and  K  emission  has  not  been  established  yet  (Shkolnik 
et  al.  2003). 

HD  187123. — This  star  may  have  (at  least)  two  planets 
orbiting  it  (Vogt  et  al.  2000).  We  provide  a  prediction  for  only 
the  inner  planet.  The  outer  planet  has  a  period  much  longer 
than  3  yr.  Consistent  with  other  planets  having  relatively  long 
periods  (d  ~  2  AU),  we  would  expect  the  flux  density  from 
this  second  planet  to  be  no  more  than  about  10~5  Jy. 

BD— 10°  3166. — This  star  has  a  planet  orbiting  it,  but  the 
distance  to  the  star  is  unknown,  although  Butler  et  al.  (2000) 
suggest  a  distance  of  less  than  200  pc.  If  that  is  the  case,  its 
flux  density  is  probably  larger  than  1 0  3  Jy. 

v  And. — There  is  uncertainty  regarding  the  mass  of  the 
outer  planet,  with  Butler  et  al.  (1999)  finding  4.6Mj,  while 
Mazeh  et  al.  (1999)  find  approximately  10Mj.  We  have  used 
the  lower  value,  realizing  that  this  factor  of  2  difference  results 
in  only  a  modest  change  in  both  vc  and  S. 
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This  star  also  shows  Ca  n  H  and  K  activity,  similar  to  that  of 
HD  179949.  However,  a  clear  link  between  the  planetary  or¬ 
bital  period  and  the  modulation  of  the  Ca  n  H  and  K  emission 
has  not  been  established  yet  (Shkolnik  et  al.  2003). 

HD  217107. — This  star  may  have  (at  least)  two  planets 
orbiting  it  (Fischer  et  al.  2001).  We  provide  a  prediction  for 
only  the  inner  planet.  The  outer  planet  has  a  period  much 
longer  than  3  yr.  Consistent  with  other  planets  having  rela¬ 
tively  long  periods  (d  ~  2  AU),  we  would  expect  the  flux 
density  from  this  second  planet  to  be  no  more  than  about 
10~5  Jy. 

HD  74156. — The  mass  of  the  outer  planet  orbiting  this 
star  is  uncertain  (in  addition  to  the  inclination  factor).  The 
cited  lower  limit  is  7.5Mj.  A  larger  value  would  lead  to  both 
a  higher  vc  and  S. 

HD  92788. — Estimates  of  the  mass  of  the  planet  orbiting 
this  star  differ  by  5%  (Udry  et  al.  2000;  Fischer  et  al.  2001). 
We  have  used  the  arithmetic  mean  of  the  masses  and  semi¬ 
major  axis  in  estimating  the  flux  density. 

HD  160691. — The  existing  radial  velocity  data  indicate,  al¬ 
though  they  do  not  provide  a  compelling  case  for,  the  presence 
of  an  outer  planet  (Jones  et  al.  2002).  If  this  second  planet  does 
exist,  its  predicted  magnetic  moment  and  radio  emission  would 
be  similar  to  those  of  HD  216435. 

e  Eli. — Quillen  &  Thorndike  (2002)  claim  the  presence  of 
an  outer  planet  based  on  the  structure  of  the  dust  disk  sur¬ 
rounding  this  star.  Even  if  such  a  planet  exists,  its  large  semi¬ 
major  axis  («40  AU)  would  imply  that  its  radio  emission 
would  be  similar  to  that  of  Neptune,  almost  certainly  at  a  fre¬ 
quency  well  below  the  Earth’s  ionospheric  cutoff. 

HD  192263. — The  existence  of  a  planet  around  this  star  was 
questioned  on  the  basis  of  photometric  variations  with  the  same 
period  as  the  radial  velocity  variations  (Henry  et  al.  2002). 
Santos  et  al.  (2003)  present  more  recent  radial  velocity  mea¬ 
surements,  which  appear  to  confirm  the  presence  of  the  planet. 

HD  47536. — The  mass  of  the  star  is  uncertain  and  in  the 
range  1-3  M0  (Setiawan  et  al.  2001).  We  have  used  the  lower 
limit  to  obtain  the  mass  of  the  planet.  Adopting  a  larger  mass 
would  imply  a  larger  predicted  flux  density,  by  a  factor  of  as 
much  as  4. 

The  star  itself  is  a  giant  (spectral  class  K1  III).  As  the 
radiometric  Bode’s  law  was  developed  for  planets  orbiting  a 
main-sequence  star,  whether  it  can  be  applied  to  planets 
orbiting  giant  stars  is  not  clear. 

HD  114762. — This  object  does  not  appear  in  Table  1,  even 
though  Farrell  et  al.  (1999)  made  a  prediction  for  its  radio 
emission,  because  it  has  since  been  reclassified  as  a  brown 
dwarf  (Schneider  2003).  Radio  emission  may  still  be  gener¬ 
ated  from  such  bodies,  though  (Berger  et  al.  2001). 

4.  VLA  74  MHz  OBSERVATIONS: 

AN  INITIAL  COMPARISON 

The  model  presented  represents  a  zeroth-order  estimate  un¬ 
der  the  assumption  that  a  solar  analog  is  valid.  To  confirm  or 
refute  the  model  some  initial  confrontation  with  observations  is 
warranted,  even  if  the  observations  cannot  be  obtained  “op¬ 
timally”  in  some  sense.  For  this  purpose  we  utilize  the  VLA 
Low-frequency  Sky  Survey  (VLSS),  an  effort  to  survey  the 
northern  sky  using  the  VLA’s  74  MHz  system  (Kassim  et  al. 
2003).  We  make  use  of  this  survey  for  two  reasons.  First,  at  a 
frequency  of  74  MHz,  the  survey  frequency  is  within  a  factor 
of  2  of  the  cutoff  frequency  of  the  radio  emission  of  Jupiter 
(~40  MHz)  and  within  a  factor  of  2  of  the  predicted  emission 
frequencies  for  nearly  one-third  of  the  current  census. 
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Fig.  2. — Field  around  70  Vir  at  74  MHz,  illustrating  the  typical  field  ob¬ 
served  in  the  VLSS.  The  noise  level  is  0.13  Jy  beam-1,  with  a  25"x25"  beam. 
The  contours  are  0.13  Jy  beam-1  times  —3,  2,  3,  4,  5,  7.07,  and  10.  The 
position  of  70  Vir  is  marked  by  a  cross.  Also  indicated  is  the  radio  source 
NVSS  J132900+140337,  likely  to  be  an  extragalactic  radio  source. 

Second,  in  producing  Table  1  we  have  reported  “burst” 
flux  densities,  which  are  a  factor  of  100  larger  than  the  median 
values.  Our  rationale  for  doing  so  is  that  we  occasionally  may 
detect  the  bursts  from  extrasolar  planets.  The  VLSS  obser¬ 
vations  are  a  series  of  overlapping  pointing  centers  on  the  sky 
with  each  pointing  center  observed  in  five  “snapshots”  of 
15  minutes  each.  While  not  as  sensitive  as  the  targeted  ob¬ 
servations  of  Bastian  et  al.  (2000),  the  VLSS  offers  the  po¬ 
tential  for  obtaining  multiple  observations  on  a  substantial 
fraction  of  the  extrasolar  planet  census.  (The  limitation  on  the 
number  of  stars  that  can  be  observed  by  the  VLSS  is  its 
southern  declination  limit  of  —30°.  Thus,  approximately  one- 
third  of  the  stars  in  the  current  census  cannot  be  observed.) 
By  observing  a  large  number  of  stars,  the  VLSS  has  the  po¬ 
tential  for  observing  a  rare,  large  burst. 

Presently,  the  VLSS  is  10%  complete,  having  surveyed  a 
0.9  sr  region.  As  an  illustration  of  its  future  utility,  as  well  as 
to  confront  the  models  with  what  limited  observations  may  be 
available,  we  have  examined  the  locations  of  the  known 
extrasolar  planets  within  the  current  survey  region.  Figure  2 
shows  the  field  around  70  Vir  as  a  typical  example. 

In  the  current  census,  five  stars  fall  within  the  region  sur¬ 
veyed  thus  far.  In  no  case  have  we  detected  emission  at  the 
location  of  an  extrasolar  planet.  Table  2  summarizes  the  upper 
limits  we  derive;  we  reproduce  the  predicted  emission  fre¬ 
quencies  and  “burst”  flux  densities  for  convenience  as  well. 
We  obtain  the  upper  limits  by  determining  the  local  noise  level 
(if),  which  is  typically  near  100  mJy.  Radio  sources  in  com¬ 
mon  to  the  74  MHz  VLSS  and  the  1400  MHz  NRAO  VLA 
Sky  Survey  (NVSS)  are  used  to  register  the  positions  of 
sources  detected  in  VLSS,  and  experience  shows  that  the  ref¬ 
erence  frames  at  the  two  frequencies  agree  to  within  approxi¬ 
mately  5".  The  radio  reference  frame  provided  by  the  NVSS 
itself  should  be  aligned  to  within  40  mas  of  the  optical  refer¬ 
ence  frame.  Thus,  the  locations  of  the  stars  hosting  extrasolar 
planets  are  known  to  a  much  better  precision  than  the  beam 
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TABLE  2 

74  MHz  VLSS  Constraints  on  Extrasolar  Planetary 
Radio  Emission 


*5model  'Lbs 


Star  (MHz)  (mjy)  (mJy) 


HD  104985  .  505  0.03  220 

HD  108874 .  54  0.06  325 

70  Vir .  546  1.6  325 

r  Boo .  60  100  275 

p  CrB .  28  16  218 


Notes. — The  model  flux  densities  are  the  “burst”  flux  densities, 
obtained  by  multiplying  the  median  levels  predicted  by  eq.  (3)  by  a 
factor  of  100.  The  upper  limits  given  by  S0bs  are  2.5  cr  in  the  VLSS  74 
MHz  images. 

(point-spread  function)  of  the  VLA  at  74  MHz  (80").  As  a 
result,  we  have  adopted  a  signal-to-noise  threshold  of  2.5  cr  in 
setting  upper  limits. 

Of  the  five  extrasolar  planets  within  the  currently  surveyed 
region  of  the  VLSS,  two  have  predicted  emission  frequencies 
well  above  the  survey  frequency  of  74  MHz.  On  the  basis  of 
the  model,  we  would  not  expect  to  see  emission  at  the  loca¬ 
tion  of  these  stars.  We  emphasize,  though,  that  one  of  our 
purposes  in  comparing  the  VLSS  observations  to  the  model 
predictions  is  an  effort  to  confront  the  model  with  whatever 
observations  are  available. 

For  the  remaining  three  extrasolar  planets,  their  emission 
frequencies  are  all  within  nearly  a  factor  of  2  of  74  MHz.  In 
these  cases  we  can  explain  the  absence  of  emission  as  due  to 
the  low  flux  densities  expected,  although  r  Boo  is  a  possible 
exception.  For  r  Boo  the  observed  upper  limit  is  within  a  factor 
of  only  3  of  the  expected  “burst”  emission  levels.  As  Farrell 
et  al.  (1999)  discussed,  r  Boo  remains  a  potentially  useful 
target  for  further  tests.  (See  also  Farrell  et  al.  2003  for  further 
observations  of  r  Boo.) 

5.  DISCUSSION  AND  CONCLUSIONS 

We  have  used  the  extrasolar  planet  census,  as  of  2003 
July  1,  and  scaling  laws  developed  in  our  solar  system  to 
predict  the  emission  frequencies  and  flux  densities  for  radio 
emission  from  extrasolar  planets.  Our  work  presents  an  update 
and  extension  of  the  initial  work  of  Farrell  et  al.  (1999). 
Figure  1  and  Table  1  present  our  primary  results.  Selection  ef¬ 
fects  due  to  the  current  planet  detection  methodology  produce 
an  apparent  deficit  of  planets  with  low  flux  densities  and  small 
emission  frequencies.  There  is  also  a  deficit  of  planets  with 
high  flux  densities  and  large  emission  frequencies  due  to  the 
well-known  deficit  of  planets  with  both  large  masses  and  small 
semimajor  axes. 

The  most  likely  planets  to  be  detected  are  those  orbiting 
nearby  stars.  Thus,  the  planet  orbiting  OGLE-TR-56  is  similar 
to  many  of  the  other  “close-in”  planets,  but  its  flux  density  is 
roughly  3  orders  of  magnitude  lower  than  most  of  the  other 
planets  in  the  census  because  of  its  much  greater  distance 
(1500  pc  vs.  a  more  typical  10  pc).  This  result  is  likely  to  be  a 
general  result  as  future  planetary  searches  move  deeper  into 
the  Galaxy. 

We  have  also  used  the  initial  observations  of  the  74  MHz 
VLSS  to  place  constraints  on  the  emission  from  a  small 
number  of  planets.  In  no  case  do  we  detect  any  emission,  and 
the  typical  upper  limit  is  approximately  0.3  Jy.  For  most  of 
the  planets,  the  predicted  flux  density  is  orders  of  magnitude 


lower  than  the  upper  limit,  but  the  upper  limit  for  the  planet 
orbiting  r  Boo  is  within  a  factor  of  2  of  the  predicted  “burst” 
flux  density. 

Most  of  the  known  extrasolar  planets  have  expected  emis¬ 
sion  frequencies  below  1000  MHz.  Indeed,  if  equation  (4)  is 
correct,  we  do  not  expect  extrasolar  planets  to  emit  strongly 
above  1000  MHz  because  their  magnetic  moments  are  not 
large  enough.  Thus,  efforts  to  detect  the  radio  emission  from 
extrasolar  planets  must  focus  on  the  low-frequency  capabili¬ 
ties  of  current  or  future  radio  telescopes. 

In  designing  an  experiment  to  detect  extrasolar  planets,  the 
potential  for  bursts  must  be  taken  into  account.  The  sensi¬ 
tivity  of  a  radio  telescope  improves  with  increasing  integra¬ 
tion  time  as  rxn~.  Longer  integration  times  improve  the 
sensitivity  but  at  the  cost  of  “diluting”  any  bursts.  For  a  burst 
of  flux  density  Sf,  and  duration  At/,,  its  average  flux  density  in 
an  integration  time  t  is  S/,(At/,/t );  i.e.,  the  burst  is  diluted  with 
time  as  t  1 .  To  the  extent  that  extrasolar  planet  radio  emis¬ 
sion  is  “bursty”  (as  is  observed  for  the  solar  system  planets), 
multiple  short  observations  are  a  better  strategy  than  long 
integrations,  as  are  employed  usually  in  radio  astronomy. 

We  now  estimate  the  sensitivity  of  various  radio  telescopes 
for  detecting  extrasolar  planets.  In  estimating  radio  telescope 
sensitivities,  particularly  for  aperture  synthesis  instruments, 
it  is  common  to  use  long  integration  times  (^8  hr).  In  contrast, 
the  magnetospheres  of  the  solar  system  planets  respond  to 
solar  wind  variations  on  timescales  of  order  1  hr.  We  shall 
use  15  minutes  as  an  integration  time  both  to  emphasize  that 
long  integration  times  may  not  be  appropriate  for  estimating 
sensitivities  and  to  allow  for  the  possibility  that  shorter 
magnetosphere  response  times  may  occur. 

The  VLA  has  been  used  already  in  attempts  to  detect  radio 
emission  from  extrasolar  planets  (Bastian  et  al.  2000;  Farrell 
et  al.  2003),  and  as  we  illustrate  here,  its  ongoing  74  MHz 
northern  sky  survey  provides  an  excellent  opportunity  to 
continue  testing  these  models.  Below  1000  MHz,  the  VLA 
operates  at  74  and  330  MHz.  Its  nominal  sensitivity  at  these 
two  frequencies  is  roughly  200  mJy/(t/15  minutes)172  and 
1.5  mJy/(t/15  minutes)172  over  bandwidths  of  1.5  and  3  MHz, 
respectively.  As  Figure  1  demonstrates,  only  if  the  radiometric 
Bode’s  law  vastly  underestimates  the  radiated  power  levels 
will  the  VLA  be  able  to  detect  radio  emissions  from  extra¬ 
solar  planets,  even  taking  into  account  the  possibility  of 
strong  bursts. 

The  Giant  Metrewave  Radio  Telescope  (GMRT)  operates 
at  frequencies  of  150,  235,  330,  and  610  MHz.  Over  this 
frequency  range,  its  nominal  sensitivity  varies  somewhat  but 
is  roughly  0.2  mJy/(t/15  minutes)172  in  a  4  MHz  bandwidth. 
Although  larger  bandwidths  (and  therefore  improved  sensi¬ 
tivities)  are  technically  possible  with  the  GMRT,  we  use  a 
4  MHz  bandwidth  in  an  attempt  to  strike  a  balance  between 
improved  sensitivity  and  avoiding  terrestrially  generated  radio 
interference.  As  Figure  1  shows,  most  extrasolar  planets  are 
below  the  detection  limit  of  the  GMRT,  unless  they  produce 
strong  bursts  (a  factor  of  roughly  100  increase  in  flux  density  ). 

The  Low  Frequency  Array  (LOFAR)  is  a  radio  telescope  in 
the  design  and  development  phase  that  is  expected  to  operate 
in  the  frequency  range  10-240  MHz.  This  frequency  range 
encompasses  the  predicted  emission  frequencies  of  roughly 
two-thirds  of  the  current  census.  The  sensitivity  of  LOFAR 
will  be  highly  frequency-  and  direction-dependent,  owing  to 
the  effect  on  the  system  temperature  from  the  Galactic  syn¬ 
chrotron  emission.  For  directions  away  from  the  Galactic 
plane,  the  design  goals  specify  a  sensitivity,  in  a  15  minute 
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integration  with  a  4  MHz  bandwidth,  around  2  mJy  at  the  lower 
frequencies  and  increasing  to  around  1  mJy  at  the  higher  fre¬ 
quencies.  There  are  a  number  of  potential  candidates  for  ob¬ 
servations  with  LOFAR.  Perhaps  most  promising  are  r  Boo, 
G1  876,  and  G1  86,  as  these  three  objects  may  be  detectable 
even  in  the  absence  of  any  flux  density  increase  due  to 
stellar  wind  variability.  Also  of  interest  are  HD  162020  and 
HD  195019,  which  may  require  only  moderate  increases  (factor 
of  10)  due  to  stellar  wind  variability  to  become  detectable. 

The  Square  Kilometer  Array  (SKA)  is  a  next-generation 
radio  telescope  that  is  expected  to  operate  above  150  MHz. 
Its  current  design  goals  suggest  that  its  sensitivity  would  be 
approximately  1  /iJy  in  a  15 -minute  integration.  This  sensi¬ 
tivity  is  sufficient  that  it  shordd  be  capable  of  detecting  the 
radio  emissions  from  the  most  massive  extrasolar  planets, 
without  relying  on  bursts  to  enhance  the  emission  levels,  and 
a  substantial  fraction  of  the  current  census  if  their  burst 
levels  are  comparable  to  what  we  find  here.  The  efficacy  of 
the  SKA  in  detecting  the  radio  emission  from  extrasolar 
planets  depends  on  two  aspects  of  its  design,  the  lowest  fre¬ 
quency  at  which  it  operates  and  its  collecting  area.  First,  if  the 
potential  SKA  frequency  range  were  to  be  increased  to  extend 
below  150  MHz,  the  number  of  planets  it  potentially  could 
detect  would  also  increase.  Second,  as  stated  currently,  the 
SKA  sensitivity  is  specified  in  terms  of  a  ratio  between  its 
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effective  collecting  area  and  system  temperature,  ^4^/7^,,. 
This  quantity  is  taken  to  be  constant.  As  Tsys  increases  below 
1000  MHz  because  of  Galactic  synchrotron  emission,  this 
design  goal  requires  that  Aeg  also  increase  below  1000  MHz 
so  as  to  keep  the  ratio  constant.  If  this  design  goal  is  relaxed 
so  that  Heff/rsys  decreases  below  1000  MHz,  the  likelihood  of 
the  SKA  detecting  the  radio  emission  from  extrasolar  planets 
also  decreases. 

All  of  our  discussions  assume  that  the  radio  emissions  from 
these  planets  would  be  driven  solely  by  the  stellar  wind  loading 
of  the  planets’  magnetospheres.  If  there  are  large,  close-in 
satellites,  like  lo  at  Jupiter,  higher  emission  frequencies  and 
flux  densities  could  result.  Also,  for  planets  that  are  close  to 
their  parent  stars,  magnetic  pressures  that  are  large  in  the  in¬ 
terior  heliosphere  may  drive  radio  emission  (Zarka  et  al.  2003). 
Future  modeling  efforts  include  integrating  the  current  census 
with  these  hot  Jupiter  emission  models. 
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